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ABSTRACT 


This  paper  presents  studies  of  the  infrared  (ir)  absorbance  and  the 
intrinsic  stress  in  thermally  grown  very  thin  films  (60A  to  700A)  of 
silicon  dioxide,  Si02-  These  data  are  combined  with  previously  obtained 

data  for  thicker  thermally  grown  films  (approximately  1300A)  to  study 
the  variation  in  intrinsic  growth  stress  close  to  the  Si/Si02  interface. 
The  combined  data  indicate  that  the  intrinsic  stress  at  Si/Si02  interfaces 
extrapolates  to  the  same  relatively  high  values  for  oxides  grown  at  700°C 
and  1000°C,  and  that  the  distribution  of  Si-O-Si  bond  angles  close  to  the 
Si/Si02  interface,  as  deduced  from  the  ir  data,  is  quantitatively  different 

than  in  the  bulk  of  the  oxide  film.  These  two  observations  are  explained  in 
terms  of  a  model  based  on  a  temperature  dependent  viscoelastic 
relaxation  of  the  oxide  stress.  This  model  emphasizes  differences  in  the 
thermal  history  of  the  Si02  near  the  Si/Si02  interface,  as  compared  to  the 
SiC>2  that  is  well-removed  from  that  interface  and  is  in  the  bulk  of  the 
film  and/or  close  to  the  'top'  surface  of  the  film.  The  observed  differences 
between  the  700°C  and  1000°C  bulk  oxides,  and  the  associated  Si/Si02 


interfaces  are  explained  in  terms  of  a  renormalized  time  scale  that  is 
defined  by  the  ratio  of  the  growth  time  to  the  viscoelastic  relaxation  time 
at  the  growth  temperature. 
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differ  by  more  than  three  orders  of  magnitude  [6],  whereas  the  times  for 
oxide  growth  of  films  of  equivalent  thicknesses  in  a  dry  ambient  at  the 
same  two  temperatures  only  differ  by  about  two  to  three  orders  of 
magnitude  [7].  These  relative  differences  in  the  two  time  scales  are  then 
shown  in  this  paper  to  promote  significant  differences  in  the  relative 
stress  relief  for  the  two  growth  temperatures  we  have  examined.  We 
explore  effects  of  this  non-uniform  thermal  history  through:  (1) 

comparisons  of  v  and  Av  as  determined  from  the  analysis  of  ir  data  in  both 
the  very  thin  and  the  thicker  oxide  films;  (2)  analysis  of  stress  and  ir  data 
in  terms  of  a  Maxwell  model  [6,8]  for  viscoelastic  relaxation  of  intrinsic 
growth  stress;  and  (3)  analysis  of  elastic  properties  of  thin  oxides  films 
using  a  combination  of  stress  and  ir  data. 

II  EXPERIMENTAL  PROCEDURES 

A.  Sample  Preparation 

The  substrate  material  used  for  oxide  growth  was  of  two  types:  2 
inch  diameter,  75  |im  thick,  2  Q-cm  p-type  silicon  wafers,  and  4  inch 

diameter,  500  |im  thick,  16-20  Q-cm  p-type  silicon  wafers.  The  75  urn 

thick  wafers  were  used  for  the  stress  measurements,  and  were  scribed 
and  broken  into  small  strips,  approximately  1  cm  by  4  cm.  The  4  inch 
diameter  wafers  were  kept  in  tact  and  subjected  to  an  HF/HNO3  based 

etch  (CP4A)  to  smooth  the  backside  of  the  wafers  and  thereby  reduce 
transmission  losses  in  the  ir  measurements  that  would  have  been  present 
due  to  the  surface  roughness  of  the  as  received  wafers.  Both  sets  of 

wafers  were  then  processed  in  the  same  way  prior  to  the  thermal 

oxidation  cycles  at  700  and  1000°C  [4].  The  wafers  were  first  subjected 
to  a  modified  RCA  clean,  and  then  were  preoxidized  to  grow  approximately 
1000A  of  oxide.  The  oxide  formed  in  this  way  was  stripped  in  buffered  HF. 
This  preliminary  oxidation  step  was  used  to  remove  contaminants  from 
the  Si  surface  prior  to  the  final  thermal  oxidations.  The  silicon  samples 
were  then  given  a  second  modified  RCA  clean,  and  were  oxidized  in  a 
double-walled,  resistiveiy-heated  fused  silica  tube  furnace  at 
temperatures  of  700° C  and  1000°C.  Silicon  wafers  were  oxidized  for 
different  times  to  produce  films  with  thicknesses  between  60A  and  700A. 
For  the  70QOC  oxides  this  corresponds  to  oxidation  times  between  about 
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20  hours  and  500  hours,  and  for  the  1000°C  oxides  to  oxidation  times 
between  about  0.1  hour  and  1.0  hours  [7].  Film  thickness  was  measured  to 
an  accuracy  of  better  than  2%  using  a  precision  manual  ellipsometer  with 
a  He-Ne  laser  (632.8  nm)  source  [1].  The  samples  used  for  the  stress 
measurements  were  then  masked,  and  the  oxide  was  stripped  from  one 
side  in  order  to  induce  curvature  from  the  stress  in  the  S1O2  film  [9].  The 

ir  measurements  were  made  on  wafers  in  which  the  thermal  oxide  was 
present  on  both  faces  [4]. 

B.  Stress  Measurements 

A  two  parallel  laser  beam  reflection  technique  was  used  to 
determine  the  total  stress  in  the  oxide  layer  as  a  function  of  oxide  film 
thickness;  details  of  this  method  have  been  discussed  in  a  paper  by  Kobeda 
and  Irene  [9].  The  technique  involves  a  measurement  of  the  radius  of 
curvature,  R,  of  the  Si  substrate  that  has  induced  the  stress  in  the  Si02 

layer.  This  curvature  is  induced  by  simply  removing  the  oxide  from  one 
side  of  the  oxidized  silicon  strip.  The  radius  of  curvature  is  then  obtained 
from  measurements  of  the  separation  between  the  two  initially  parallel 
laser  beams  reflected  from  the  oxidized  surface  of  the  Si  sample. 
Calibration  is  achieved  by  comparisions  with  flat  mirror  surfaces,  and 
with  mirrors  of  known  curvature,  R  -  20  m  and  50  m.  These  two  mirror 
radii  are  within  the  range  of  curvatures  present  in  most  of  the  stressed 
Si/Si02  heterostructures.  The  total  stress,  of,  in  the  Si02  film  is 

determined  from  a  standard  beam  bending  relationship  [9,10]; 

Of  -  (E'  ts2)/(6(1-/)tfR);  (1) 

where  ts  and  tf  are  the  substrate  and  film  thicknesses,  respectively  and  E' 
and  y'  are  Young's  modulus  and  Poisson's  ratio  for  Si,  respectively.  Since 
the  thermal  expansion  coefficients  of  Si  and  Si02  differ  significantly  [9], 
there  is  an  appreciable  contribution  to  of  from  thermal  effects,  i.e., 
stress  that  is  induced  by  virtue  of  film  growth  at  the  elevated 
temperatures,  700°C  and  1000°C,  and  stress  measurements  made  at  room 
temperature,  approximately  20°C.  The  total  film  stress  is  of  »  oth  +  <*j. 
where  oth  is  a  thermally  induced  stress  that  can  be  calculated  from  the 
known  values  of  the  thermal  expansion  coefficients  of  Si  and  Si02,  and  oj 
the  intrinsic  growth  stress,  i.e.,  the  stress  that  is  generated  through  the 
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film  growth  process.  This  intrinsic  growth  stress  is  the  focal  point  of 
this  study. 

We  have  been  able  to  extend  stress  measurements  to  oxide 
thicknesses  in  the  range  of  100A  by  the  use  of  very  thin  75  |im  Si 
substrates.  Since  Eqn.  (1)  includes  scaling  factors  of  ts2  and  1/R,  a 
thinner  substrate  will  yield  a  smaller  radius  of  curvature,  i.e.,  more 
curvature  for  a  given  stress  level  and  film  thickness.  The  reported  values 

of  stress  in  this  paper  are  generally  the  average  of  four  measurements 

made  on  each  of  three  different  samples.  This  multiplicity  of  stress 
measurements  has  enabled  us  to  determine  oj  to  an  average  uncertainty  of 

about  0.5  x  109  dynes/cm2. 

C.  Infrared  Spectroscopy  Measurements 

The  ir  absorbance  measurements  reported  here  were  obtained  as 
described  in  Refs.  1-5.  The  measurements  on  the  very  thin  oxides  involve 

averages  over  several  different  samples.  We  have  found  that  in  some 

instances,  noise  reduction  in  the  data  through  multiple  point  smoothing 
can  lead  to  distortions  of  the  spectra  that  become  evident  in  samples  such 
as  these  which  have  asymmetric  absorption  features.  All  of  the  data  in 
this  paper,  except  for  the  60A  samples,  rely  on  noise  reduction  through 
long  term  data  acquisition  rather  than  any  form  of  post  acquision 
processing.  We  have  taken  care  to  eliminate  distortion  in  the  60A  samples 
by  the  judicious  choice  and  application  of  smoothing  routines.  Figure  1 
indicates  the  ir  absorbance  for  the  bond-stretching  vibration  at  about 
1075cm*1  for  a  thermally  grown  Si02  thin  film.  There  are  two  other 
first-order  ir  active  vibrations  in  Si02,  a  bond-bending  band  at  about 
800cm'1,  and  a  bond-rocking  band  at  about  460cm*1  [1].  We  have 
previously  reported  systematic  changes  in  the  1075cm*1  band  with 
different  growth  and  annealing  temperatures  [1-5].  Changes  in  the 
character  of  other  two  bands  are  significantly  smaller  and  therefore  more 
difficult  to  use  to  monitor  the  effects  of  different  growth  and/or 
annealing  temperatures.  The  insert  in  the  figure  shows  the  relative  motion 
of  the  oxygen  atom  for  the  bond-stretching  vibration;  this  motion  is  in  the 
plane  of  the  Si-O-Si  bond,  and  in  a  direction  parallel  to  a  line  joining  the 
two  silicon  atoms. 
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In  this  report,  we  restrict  our  observations  and  the  discussions 
following  to  the  center  frequency,  v,  and  the  full-width  at  half-maximum 
(FWHM),  Av,  of  the  bond-stretching  feature;  these  are  also  indicated  in  Fig. 
1.  The  values  for  v  and  Av,  reported  in  the  ramainder  of  this  paper, 
represent  an  average  over  six  spectra  taken  on  samples  from  the  same 
furnace  run.  The  experimental  uncertainties  in  v  and  Av  come  primarily 
from  two  sources,  the  signal  to  noise  ratio  in  the  data  acquisition 
process,  and  the  inherent  resolution  of  the  ir  spectrophotometer  for  the 
particular  set  of  conditions  under  which  the  data  was  taken.  The  signal  to 
noise  ratio  is  the  dominant  factor  in  the  new  data  that  we  are  reporting 
for  the  thin  thermally  grown  oxides,  and  also  varies  with  film  thickness. 
For  the  thinnest  films,  60A,  the  uncertainties  are  approximately  1  cm'1 
for  both  v  and  Av,  where  as  for  thicker  films,  greater  than  about  500A,  the 
respective  uncertainties  are  0.6cm-1  for  v,  and  0.2cm-1  for  Av  [4]. 

Ill  EXPERIMENTAL  RESULTS 

We  have  obtained  new  data  for  v,  Av  and  oj  for  thin  oxide  films  grown 

at  either  700°C  or  1000OC.  In  Figs.  2-6,  we  have  combined  these  new  data 
with  data  previously  reported  [2-5]  for  thicker  thermally  grown  Si02 

films  in  order  to  emphasize  the  new  features  that  emerge  as  we  probe  the 
oxide  properties  closer  to  the  Si/SiC>2  growth  interface. 

Figures  2-6  then  present  the  results  of  ir  (v  and  Av)  and  stress  (o j) 
measurements  for  films  with  thicknesses  between  about  60A  and  1300A. 
The  data  shown  in  Figs.  2a  and  b  were  taken  on  thick  oxides 
(approximately  1300A  thick)  and  indicate  the  way  that  v  and  Av  vary  with 
the  growth  temperature,  Tox.  Specifically,  v  increases  as  Tox  is 
increased,  and  Av  decreases  as  T0x  is  increased.  The  results  shown  in  Fig. 
3  for  Av  as  a  function  of  v  were  also  taken  on  relatively  thick  oxides, 
average  thickness  1 300A,  and  indicate  a  correlation  between  Av  and  v  that 
is  obeyed  universally  by  all  of  these  thick  oxides  that  we  have  studied  [2- 
5].  These  thin  film  layers  include:  (1)  thermal  oxides  grown  in  dry  and 
steam  ambients  at  temperatures  between  800°C  and  1150°C;  (2)  thermal 
oxides  that  have  been  subsequently  annealed  at  temperatures  up  to 
1200°C;  (3)  oxides  deposited  by  remote  plasma-enhanced  chemical-vapor 
deposition  (RPECVD)  between  100°C  and  550°C;  (4)  and  oxides  deposited 
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by  RPECVD  and  then  annealed  at  temperatures  up  to  1000°C.  The  solid  line 
is  a  linear  regression  analysis  and  is  shown  to  quantify  the  correlated 
behavior  of  Av  and  v  for  these  thick  oxides. 

Figure  4  presents  the  variation  of  the  center  frequency,  v,  of  the  Si- 
O  bond-stretching  vibration  as  a  function  of  film  thickness  for  oxides 
grown  at  700°C  and  1000°C.  The  samples  grown  at1000°C  have  larger 
values  of  v  at  any  given  film  thickness  than  the  samples  grown  at  700° C. 
In  all  cases,  however,  as  the  film  thickness  is  decreased,  v  also  decreases. 
Moreover,  as  the  film  thickness  goes  to  zero  both  sets  of  data  extrapolate 
to  approximately  the  same  frequency  for  the  bond-stretching  vibration, 
1055cm"1±  Icnrrl.  The  solid  lines  in  the  figure  are  derived  from  a 
logarithmic  regression  analysis.  Figure  5  displays  the  corresponding 
behavior  for  the  variation  of  Av  as  a  function  of  film  thickness  for  oxides 
grown  at  700°C  and  1000°C.  For  the  oxides  grown  at  700°C,  the  value  of 
Av  is  relatively  constant  at  about  79  cm'1,  while  for  the  oxides  grown  at 

1000°C,  Av  decreases  significantly  in  the  thinnest  oxides.  The  data  of  Fig. 

5  then  demonstrate  a  break-down  in  the  universal  relationship  between  v 
and  Av  in  the  thinnest  oxide  films  that  we  have  studied. 

Figure  6  includes  a  plot  of  intrinsic  growth  stress  versus  film 
thickness  for  oxides  grown  in  a  dry  ambient  at  700°C  and  1000°C  and  for 
film  thicknesses  to  approximately  100A  [6].  At  any  given  film  thickness, 
the  oxides  grown  at  700° C  have  a  higher  intrinsic  growth  stress  than 
those  grown  at  1000°C,  with  the  stress  level  in  both  sets  of  oxides 

increasing  as  the  oxide  thickness  is  decreased.  The  data  at  each  growth 

temperature  converge  to  the  same  value  of  the  intrinsic  stress, 
approximately  4.6  x  109  dynes/cm2,  as  the  oxide  thickness  goes  to  zero. 
This  occurrence  of  a  common  maximum  value  of  intrinsic  growth  stress  at 
the  Si/Si02  interface  is  not  surprising  since  the  origin  of  this  interfacial 

stress  is  the  constant  molar  volume  change  strain  that  occurs  at  the  Si 
surface  during  the  oxidation  process  [6].  It  should  also  be  noted  from  Fig. 
4,  that  in  addition  to  the  extrapolated  values  of  <rj  being  the  same  for 
films  grown  at  700°C  and  1000°C,  the  extrapolated  values  of  v  are  also 
the  same.  We  will  consider  the  implications  of  this  in  our  discussion  of 
the  experimental  data  that  follows. 
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IV  DISCUSSION 

We  now  proceed  to  discuss  the  data  presented  above  in  the  context 
of  a  model  for  the  thermal  oxidation  of  silicon  that  emphasizes  an  aspect 
of  that  growth  process  that  we  characterize  as  a  non-homogeneous 
thermal  history.  The  model  is  based  on  a  comparison  between  two  time 
scales:  (1)  the  time  required  to  grow  an  oxide  film  of  a  specified 
thickness;  and  (2)  the  viscoelastic  relaxation  time  at  that  same  growth 
temperature.  The  notion  of  viscoelastic  relaxation  arises  from  the  fact 
that  there  is  a  considerable  increase  in  molar  volume  between  the 
underlying  silicon  material  and  the  growing  oxide  [6];  this  in  turn  is 
responsible  for  the  build  up  of  a  compressive  stress  in  the  oxide  film  in 
the  immediate  vicinity  of  the  Si/Si02  interface.  The  stress  at  a  given 

point  in  the  oxide  film  can  be  relieved  during  the  remainder  of  the  oxide 
growth  process  only  if  the  time  for  viscoelastic  relaxation  at  the  growth 
temperature  is  much  shorter  than  the  time  required  for  completion  of  the 
oxidation  process.  There  are  two  limiting  cases  that  we  wish  to  compare 
in  a  qualitative  manner. 

The  first  case  is  concerned  with  oxide  growth  at  a  relatively  low 
temperature  at  which  the  viscoelastic  relaxation  time  is  longer  than  the 
oxidation  time.  Under  these  conditions,  there  is  very  little  viscoelastic 
relaxation  and  any  stress  generated  at  the  growth  surface  is  not  relaxed 
effectively  during  the  subsequent  film  growth.  Under  these  conditions 
there  will  be  a  relatively  small  stress  gradient  between  the  growth 
interface  and  the  surface  of  the  oxide  film;  i.e.,  between  the  "newest"  and 
"oldest"  regions  of  the  oxide  film.  The  other  extreme  is  at  a  higher 
oxidation  temperature  where  the  time  for  viscoelastic  relaxation  can  be 
comparable  to,  or  even  shorter  than  the  oxidation  time.  Under  these 
conditions  the  "new"  oxide  generated  at  the  Si/SiC>2  interface  is  under 

stress  due  to  the  change  in  molar  volume,  but  "older"  oxide  in  the  rest  of 
the  oxide  layer  can  experience  a  substantial  degree  of  stress  relief  via 
viscoelastic  relaxation.  We  now  demonstrate  that  the  two  cases  discussed 
above  are  realized  in  oxides  grown  in  dry  oxygen  ambients  at  temperatures 
of  700°C  and  1000°C,  respectively. 
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In  order  to  study  the  distribution  of  intrinsic  stress,  and  the 
associated  strain  in  thermally  grown  oxide  films,  we  determined  the 
stress  as  a  function  of  film  thickness  in  oxides  grown  at  different 
temperatures  that  span  a  significant  range  of  growth  rate  and 
viscoelastic  relaxation  rate  time  scales,  in  addition,  we  have  also  studied 
the  center  frequency,  v,  and  FWHM,  Av,  of  the  Si-0  bond-stretching  mode 
and  have  established  in  previous  publications  that  the  center  frequency  of 
this  mode  gives  direct  information  about  the  bond-angle,  20,  at  the  oxygen 
atom  sites  [1-5],  and  that  relief  of  intrinsic  stress  is  accompanied  by 
systematic  increases  in  v,  and  therefore  in  20,  and  also  by  decreases  in  Av 
[1-5]. 


The  systematic  variations  in  the  stretching  frequency,  v,  and  the 
intrinsic  growth  stress,  <jj,  shown  respectively  in  Figs.  4  and  6  for  the 
films  grown  at  1000°C,  demonstrate  a  continuous  variation  in  the  thermal 
history  from  highly  stressed  oxide  at  the  growth  interface,  as  displayed  in 
the  thinnest  films,  to  a  more  relaxed  and  essentially  annealed  oxide 
characteristic  of  the  bulk,  as  displayed  in  the  thicker  films.  Since  the 
stress  in  the  Si02  film  near  the  growth  interface  is  in  the  primarily  in  the 

the  the  plane  defined  by  that  interface  and  defined  here  as  the  x-y  plane, 
the  effective  stress  in  the  z-direction,  or  the  direction  of  film  growth, 
must  be  relatively  tensile.  This  follows  from  the  tendency  of  a  material 
to  retain  a  constant  volume,  and  is  the  concept  upon  which  Poisson's  ratio 
is  defined.  According  to  the  simple  Maxwell  model  of  a  solid  [6],  we  expect 
stress  relaxation  to  occur  by  viscous  motion  of  the  Si02  material  away 
from  the  growth  interface.  The  out-of-plane  strain  relaxation  rate,  de/dt, 
in  the  Maxwell  model  is  then  proportional  to  the  in-plane  stress: 

de/dt  »  oj/ri;  (2) 

where  n  is  the  oxide  viscosity.  The  smaller  the  value  of  the  viscocity  the 
higher  the  strain  relaxation  rate,  and  consequently  the  more  rapid  the 
relaxation  of  the  intrinsic  growth  stress  in  the  region  removed  from  the 
Si/S i02  interface.  The  equilibrium  viscosity  is  an  exponential  function  of 

temperature,  T,  as  shown  in  the  following  equation  [8]; 


h  -  ho  ®xP(Erj/RgT); 


(3) 
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where  is  a  constant  representing  the  minimum  viscocity,  is  an 
empirically  determined  activation  energy  and  Rg  is  the  gas  constant.  If  the 
rate  of  stress  relief  is  assumed  to  be  proportional  to  the  stress,  as  is 
usually  the  case,  then; 

doj/dt  ■  -  (1/x)  oj.;  (4) 

where  t  is  a  constant,  namely  the  viscoelastic  relaxation  time.  This 
relaxation  time  is  assumed  to  be  independent  of  the  time  scale  of  the 
experiment.  One  can  then  integrate  Eqn.  (4)  and  obtain  the  following 
expression  for  the  stress  relief  process; 

cj(t)  -  <rj(o)  exp  (-  t/x);  (5) 

where  oj(o)  is  the  maximum  intrinsic  stress  at  the  growth  interface,  and  t 

is  the  oxidation  or  annealing  time.  From  the  experimental  data  displayed 
in  Fig.  6,  we  find  an  empirically  determined  value  for  oj(o)  of 

approximately  4.6  x  109  dynes/cm2.  Relaxation  times  have  been  predicted 
from  viscoscity  measurements  by  assuming  that  the  relaxation  time,  x,  is 
proportional  to  the  viscosity  [6]; 

x-r|/G,  (6) 

where  G  is  the  shear  modulus.  The  shear  modulus  is  a  weak  function  of 
temperature,  except  for  the  significant  decrease  in  value  near  950- 
1000°C  [11],  so  that  the  very  strong  exponential  temperature  dependence 
of  n  dominates  over  most  of  the  temperature  range  [6].  In  order  to 
estimate  relaxation  times  from  the  data  displayed  in  Fig.  6,  it  is 
necessary  to  take  into  account  that  the  measured  values  of  stress 
represent  an  average  integrated  over  the  entire  film  thickness.  Therefore, 
in  order  to  be  consistent  with  this  aspect  of  the  experimental  data,  the 
value  of  t  used  in  Eqn.  (5)  and  obtained  from  experimental  data,  must  be  an 
average  oxidation  time  integrated  over  the  film  thickness.  Such  an  average 
value  for  t  can  be  obtained  by  integrating  an  empirical  relationship  for 
thickness  as  a  function  of  oxdiation  time,  e.g.,  the  Grove-Deal  model  [12], 
and  normalizing  this  integral  through  division  by  the  total  film  thickness. 
Relaxation  times  estimated  from  this  procedure  appear  in  Table  I.  It  is 
also  possible  to  estimate  relaxation  times  from  the  ir  data,  where  we 


have  already  shown  that  the  center  frequency  of  the  bond-stretching 
vibration  provides  a  measure  of  the  bond-angle  at  the  oxygen  atoms  sites. 
We  develop  this  argument  below,  and  indicate  the  validity  of  an 
assumption  that  is  based  on  stress  relaxation  through  changes  in  that 
bond-angle  [5]. 

The  systematic  relationships  oetween  v  and  Av  displayed  in  Figs.  2 
and  3  have  been  explained  in  previous  publications  in  terms  of  a  model 
based  on  the  local  atomic  structure,  in  particular  the  bond-angle,  20,  at 
the  bridging  oxygen  atom  sites  (see  Fig.  1)  [1,2].  Within  the  framework  of 
this  model,  the  average  or  central  frequency  of  the  ir  active  bond¬ 
stretching  vibration  is  given  by: 

v  «■  v0  sin  9,  where  v0  =  1116.5  cm-1;  for  0-75  degrees  (7) 

This  relationship  between  v  and  0  is  based  on  a  central  force  model  for  the 

vibrational  modes  of  Si02  [13].  The  model  equation  neglects  motion  of  the 

oxygen  atoms,  and  considers  only  two-body  central  bond-stretching 
forces,  neglecting  specifically  the  smaller  three  body  bond-bending 
forces,  etc.  The  neglect  of  these  terms  renormalizes  the  force  field  with 
the  consequence  of  defining  v0  empirically,  rather  than  calculating  it  from 

microscopic  forces,  mode  masses,  etc.  In  any  case  the  functional 
dependence  of  v  on  0  is  the  important  aspect  of  this  model,  and  this  is  not 
compromised  by  the  approximation  inherent  in  Eqn.  (7).  The  central  force 
model  also  neglects  small  changes  in  the  Si-0  bond-length  with  changes  in 
the  bond-angle  [14].  These  changes  in  the  bond-length  are  implicitly 
included  in  the  empirical  parameter,  v0.  It  should  be  noted  that  similar 

force  constant  approximations  based  on  central  forces  alone,  do  not  apply 
for  the  other  two  first  order  vibrational  modes,  so  that  any  systematic 
changes  in  0  with  growth  and  or  annealing  temperatures  are  not  expected 
to  have  similar  systematic  effects  on  the  frequencies  of  these  other  two 
first-order  modes.  This  is  consistent  with  our  observation  that  changes  in 
the  bond-bending  and  rocking  modes  with  growth  and/or  annealing 
conditions  are  smaller,  and,  as  such  it  has  not  proven  possible  to  use  them 
as  parameters  for  tracking  relief  of  intrinsic  growth  stress.  We  use  X-ray 
data,  and  assign  the  average  Si-O-Si  bond-angle  a  value  of  150  degrees  for 
an  oxide  grown  at  1150°C.  This  choice  of  bond-angle  is  not  critical  in  the 


12 


development  that  follows  in  the  sense  that  it  represents  a  bench  mark 
position  which  simply  serves  as  a  point  of  reference  for  the  comparisons 
that  are  made.  The  choice  of  1150°C  represents  a  temperature  at  which 
viscoelastic  relaxation  time  is  short  compared  to  the  time  scale  for 
growth.  For  a  relatively  thick  oxide  (1300A)  grown  at  1150°C  we  find  a 
vibrational  frequency  of  1078.5cm-1.  Combining  this  with  a  bond-angle  of 
150  degrees  yields  the  value  of  v0  given  in  Eqn.  (7). 

At  this  point,  we  note  that  the  average  Si-Si  (2nd  neighbor) 
distance,  dsi-Si  *s  a*so  a  function  of  0,  through  the  relation: 

dsi-Si  -  2  ro  sine  (8) 

where  r0  is  the  Si-O  bond-length;  r0  has  been  shown  to  be  a  slowly  varying 
function  of  0  [14],  with  the  result  that  the  assumption  of  a  constant  r0  is 
consistent  with  the  approximations  used  in  Eqn.  (7).  The  representation  of 
the  Si-Si  distance  through  Eqn.  (8)  is  important  in  our  discussion  of  stress 
relief,  since  it  has  been  shown  that  the  density  of  Si02  thermally  grown 
oxide  films  scales  inversely  with  [dsi-Si]3  [1-5].  This  has  been 

established  through  correlated  trends  in  the  behavior  of  the  index  of 

refraction,  n,  of  SiC>2  (in  a  transparent  region  of  the  spectrum,  632.8nm) 
with  the  vibrational  frequency  v.  Since  v  is  a  function  of  0  through  Eqn.  (7), 
and  since  n  is  proportional  to  density  [1],  a  correlation  between  v  and  n  is 
expected,  i.e.,  n  increases  as  v  decreases.  Using  Eqn.  (7)  and  (8),  we  have 
been  able  to  quantify  that  correlation  by  predicting  the  scale  factor  for 
the  approximately  linear  relationship  observed  between  v  and  n  [1].  The 
implication  of  a  larger  bond-stretching  frequency  is  a  larger  average 
bond-angle  at  the  oxygen  atom  sites,  and  hence  a  larger  Si-Si  distance  and 
a  reduced  density.  As  will  be  discussed  later,  the  changes  in  the  Si-O-Si 
bond  angle  and  density  are  themselves  a  manifestation  of  stress  relief.  We 

have  developed  an  additional  aspect  of  our  model  by  interpreting  the 

correlated  data  in  Fig.  3  for  bulk  oxides  and  then,  in  Fig.  7,  superimposing 
data  from  the  thin  oxides  to  highlight  how  the  local  atomic  structures 
changes  as  the  Si/Si02  interface  is  approached. 

We  have  shown  that  the  width  of  the  Si-0  stretching  absorption  band 
in  the  thick  oxides  could  be  understood  in  terms  of  a  statistical 
distribution  of  Si-O-Si  bond-angles.  If  Av  comes  from  such  a  summation 
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over  narrow  vibrational  bands,  each  associated  with  a  different  Si-O-Si 
bond  angle,  then  Av  -  dv  and  we  can  simply  differentiate  Eqn.  (7)  to  obtain 
the  functional  relationship  between  Av  and  9: 

dv«v0cos9de.  (9) 

Dividing  Eqn.  (9)  by  Eqn.  (7),  and  then  dividing  and  muliplying  the  right- 
hand  side  of  that  relationship  by  0  yields  the  following: 

dv/v  *  0  cot0  (d0/0)  (10) 

We  have  replotted  the  data  of  Fig.  7  in  Fig.  8  using  the  representation  of 
Eqn.  (10).  The  slope  derived  from  this  replotting  of  the  Fig.  3  data  is  dO/0 - 
0.22;  this  value  of  the  average  bond-angle  spread  is  in  excellent 
agreement  with  X-ray  data  [15].  We  have  also  included  data  from  the 
thinnest  oxides,  where  we  find  that  the  slope,  d0/0,  increases  as  the  oxide 
thickness  is  reduced.  This  will  be  discussed  in  more  detail  later  on  in  the 
paper.  We  first,  however,  return  to  the  data  in  Fig.  4  where  we  have 
displayed  v  as  a  function  of  film  thickness  for  oxides  grown  at  700° C  and 
1000°C. 

In  the  context  of  the  central  force  model  for  v,  decreases  in  v  with 
decreasing  film  thickness  also  reflect  decreases  in  9.  From  Eqn.  (8) 
decreses  in  9  also  imply  a  decreasing  Si-Si  distance  and  hence  a  more 
dense  or  more  compacted  film.  This  is  consistent  with  the  data  that  is 
displayed  in  Fig.  6  which  shows  increases  in  intrinsic  compressive  stress 
as  the  film  thickness  is  reduced.  The  data  in  Fig.  6  are  particularly 
interesting  with  regard  to  the  extrapolated  values  of  stress  for  zero  film 
thickness.  The  observation  that  the  data  for  oxides  grown  at  both  700°C 
and  1000°C  extrapolate  to  essentially  the  same  value,  then  leads  us  to 
propose  that  the  constant  value  of  stress  is  primarily  a  manifestation  of 
the  increase  in  molar  volume  that  occurs  in  the  conversion  of  Si  to  Si02 
[6]. 


The  data  in  Fig.  5  indicates  that  Av  remains  relatively  constant  as  a 
function  of  film  thickness  for  the  oxides  grown  at  700° C,  but  that  Av 
decreases  dramatically  for  the  very  thinnest  oxides  grown  at  1000°C. 
These  observations,  coupled  with  the  data  displayed  in  Fig.  4  indicate  that 
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as  oxide  thickness  is  reduced  the  linear  relationship  between  Av  and  v, 
displayed  in  Figs.  3  for  thick  oxides,  starts  to  break  down.  This  is  evident 
from  the  data  displayed  in  Fig.  7,  where  we  have  plotted  Av  versus  v  for 
the  thin  oxides.  We  have  also  included  in  Fig.  7  the  result  of  a  linear 
regression  analysis  of  the  thick  oxide  data  of  Fig.  3.  The  data  for  the  thin 
oxides  fall  below  the  average  behavior  of  the  thick  oxides  indicating  that 
the  center  frequency  of  the  absorption  band  has  shifted  to  lower 
wavenumbers  relative  to  the  distribution  of  bond-angles.  If  the  data  from 
the  thin  oxides  had  the  same  slope  as  the  regression  analysis  line  but  was 
displaced  by  a  constant  amount,  then  we  would  expect  smaller  values  of  v 
to  have  correspondingly  smaller  values  of  Av.  But,  since  the  relative 
distribution  of  bond-angles,  AO/0,  depends  both  on  the  spread  in  bond- 
angles,  A0,  and  the  value  of  0  itself  [see  Eqn.  (10)],  a  slightly  smaller 
value  of  Av  does  not  necessary  imply  a  smaller  relative  distribution  of 
bonds  angles.  In  fact  we  shall  demonstrate  just  the  opposite  exists.  But 
before  proceeding,  it  should  be  noted  that  Boyd  [16]  has  also  reported 
similar  trends  in  the  ir  absorption  spectra  of  thin  oxides,  but  his 
observations  were  based  on  direct  transmission  data,  rather  than 
calculated  values  of  the  absorbance,  and,  as  such  any  quantitative 
comparisons  with  the  trends  he  has  reported  are  not  possible.  Another  way 
of  looking  at  the  behavior  displayed  in  Fig.  8  is  based  on  the  changes  in  the 
local  atomic  structure,  specifically  the  Si-O-Si  bond  angle  20,  as  the 
growth  interface  is  approached.  This  can  be  treated  in  a  quantitative 
manner  by  including  the  thin  oxide  data  in  a  plot  of  dv/v  versus  0  cot0,  and 
then  using  the  slope  of  that  plot  as  a  basis  for  determining  a  value  for 
d0/0.  This  is  shown  in  Fig.  8,  where  we  have  included  the  regression 
analysis  fit  to  the  thick  oxide  data  as  a  reference,  and  then  displayed 
representative  thin  oxide  data.  We  find  a  value  of  d0/0  -  0.47  for  the 
thinnest  oxides.  This  is  significantly  greater  than  the  value  of  d0/0  -  0.22 
for  the  thick  oxides,  and  in  fact  is  too  large  to  have  any  direct  physical 
significance;  i.e.,  the  model  used  to  analyze  the  data  no  longer  applies. 
Specifically,  we  interpret  this  large  value  of  d8/0  to  mean  that  the 
distribution  of  bond-angles  near  the  Si/Si02  interface  is  no  longer 
statistical  in  nature,  and  is  determined  instead  by  other  types  of  bonding 
constraints  associated  with  the  differences  in  molar  volume  between 
crystalline  silicon  and  Si02<  and  therefore  indicates  that  the  spread  in  the 
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bond-angle  at  the  oxygen  atom  sites  is  being  increased  by  the  build-up  of 
growth  stress  as  the  Si/Si02  interface  is  approached.  At  first  glance  a 

seemingly  greater  relative  spread  in  bond-angles  appears  to  be  counter¬ 
intuitive,  since  one  would  have  expected  constraints  in  the  local  atomic 
structure,  developed  from  the  compressive  stress  at  the  growth  interface, 
to  always  decrease  the  spread  in  Si-O-Si  bond  angles.  This  point  clearly 
needs  additional  study,  both  experimentally  and  through  modelling  of  the 
way  the  intrinsic  stress,  generated  by  the  strain  at  the  growth  interface, 
is  relieved  during  subsequent  growth  and  annealing. 

The  ir  data  obtained  for  the  bond-stretching  vibrations  and  displayed 
in  Fig.  4  can  also  be  analyzed  in  terms  of  a  Maxwell  model  for  viscoelastic 
relaxation.  We  assume  a  relationship  of  the  form: 

vmax  ■  v  *  (vmax  *  vmin)  exP  (“*/*)  (1 1 ) 

where  vmax  is  determined  from  the  high  temperature  (1150°C)  thick 

(1300A)  oxide  data  and  is  set  equal  to  1079cm*1 ,  the  approximate 
asymptote  in  the  v  versus  growth  temperature  plot  [see  Fig.  2];  vmjn  is 

estimated  from  the  thinnest  oxides  (60A)  grown  at  700°C,  and  has  a  value 
of  approximately  1055cm*1.  As  in  Eqn.  (5),  t  is  an  average  growth  (or 
annealing)  time  and  x  is  a  viscoelastic  relaxation  time.  We  have  calculated 
values  of  x  from  the  ir  data  and  these  are  displayed  in  Table  I  along  with 
those  obtained  from  the  stress  data.  The  relaxation  times  computed  for  a 
growth  temperature  of  1000°C,  based  on  the  data  from  a  607A  oxide  film 
are  in  very  close  agreement  with  an  average  value  for  x  of  approximately 
0.75  hours  or  equivalently,  approximately  2.7  x  10^  s.  The  data  for  700°C 
growth,  based  on  comparision  with  an  oxide  film  that  is  701 A  show  a 
somewhat  greater  spread  (about  ±  50%)  with  an  average  value  of  about 
800  hours  or  approximately  2.9  x  10^  s.  Values  of  x  estimated  from 
viscocity  and  shear  modulus  data  [6]  are  also  displayed  in  Table  I.  There  is 
order  of  magnitude  agreement  between  the  two  methods  of  estimating  x  at 
1000°C,  0.2  hours  from  the  viscosity/shear  modulus  approach  as  compared 
to  0.65  hours  from  Eqn.  (11).  However,  the  viscosity  data  yields  a 
significantly  higher  value  of  x  at  700°  C,  approximately  1270  hours  as 
compared  to  332  hours  from  Eqn.  (11)  [6].  We  suspect  the  differences 
between  these  determinations  for  x  at  700°C  arise  at  least  in  part  from 
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possible  errors  in  the  extrapolation  of  viscosity  and  shear  modulus  data. 
We  now  compare  the  relative  time  scales  for  relaxation  and  growth,  and 
discuss  our  observations  relative  to  the  different  stress  gradients  in  the 
films  grown  at  700°C  and  1000OC  in  terms  of  these  two  time  scales. 

We  note  the  times  required  to  grow  thermal  oxides  in  a  dry  oxygen 
ambient  at  1000°C  range  between  0.25  and  2.0  hours  for  oxide  thickness 
between  200A  and  lOOOA,  whilst  times  between  100  and  500  hours  are 
required  for  the  same  thickness  range  at  700° C.  For  example,  the  growth 
time  was  approximately  1.2  hour  for  a  700A  oxide  grown  at  1000OC  , 
whereas  for  a  600A  oxide  grown  at  700° C,  the  time  was  approximately 
350  hours.  This  means  that  the  time  scale  for  the  growth  of  the  oxide  at 
1000°C  would  allow  relaxation  of  the  "old  oxide”  in  thicker  films,  but  not 
the  "newer  oxide”  near  the  growth  interface  [the  viscoelastic  relaxation 
time  at  1000°C  is  approximately  0.75  hours].  In  contrast  in  the  films 
grown  at  700°  C,  the  growth  time  scale  is  too  short  for  appreciable 
relaxation  [the  average  viscoelastic  relaxation  time  at  700°C  is  about  800 
hours].  These  differences  in  the  relative  growth  and  relaxation  times  are 
clearly  visible  in  the  data  displayed  in  Figs.  4  and  6.  Implicit  in  the 
estimates  of  x  that  we  have  made  from  the  stress  and  ir  data  is  an 

assumption  that  the  ir  data  reflects  in  some  way  the  strain  that  is 

associated  with  the  growth  induced  intrinsic  stress.  This  is  assumption  is 
supported  through  an  analysis  that  is  given  below. 

One  way  to  demonstrate  that  the  ir  data  and  stress  data  are 
complementary  in  giving  information  about  stress  and  stress  relief  is  to 
demonstrate  that  the  ir  data  and  stress  data  can  be  combined  to  yield  one 
of  the  elastic  constants  of  Si02-  We  have  previously  demonstrated  a 

proportionality  between  the  average  intrinsic  stress  in  thick  oxide  films 
and  the  strain  as  estimated  from  the  ir  data  [4].  We  now  present  the 
assumptions  used  in  that  analysis  and  extend  it  to  include  the  new  thin 
oxide  data.  Using  the  model  described  in  Eqns.  (7)  and  (8),  we  proceed  in 

the  following  manner.  We  first  assume  that  changes  in  the  Si-Si  distance 

are  an  appropriate  relative  strain  parameter  and  that  these  may  be 
deduced  from  the  ir  data  via  the  dependence  of  both  v  and  dsi-Si  on  We 
therefore  define  the  relative  strain,  e,  as  given  by: 
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e  -  Adsi-Si  '  ^Si-Si-  (12) 

We  define  20i  as  the  Si-O-Si  bond-angle  in  an  oxide  grown  at  1150°C,  and 
and  r0  as  the  corresponding  Si-0  bond-length.  If  we  neglect  any  changes  in 
the  Si-O  bond-length  with  changes  in  the  Si-O-Si  bond  angle,  then  if  202 
as  the  Si-O-Si  bond-angle  in  an  oxide  grown  at  a  lower  temperature;  e.g., 
700°C.  Substitution  of  Eqn.  (8)  into  Eqn.  (12)  yields  the  following 
approximate  value  for  the  strain: 

e  -  (2r0  sine-j  -  2r0  sin  02)/2ro  sin  0-j.  (13) 

We  have  previously  shown  that  it  was  valid  to  neglect  any  changes  in  the 
Si-O  bond-length  when  considering  the  applicability  of  the  central  force 
model  for  relating  the  Si-0  bond-stretching  frequency  to  the  Si-O-Si 
bond-angle.  However,  for  the  purposes  of  estimating  the  strain,  these 
small  changes  in  the  Si-0  distance  must  be  taken  into  account.  Consider 
first  the  approximation  given  by  Eqn.  (13),  where  the  changes  in  the  Si-0 
bond-length  are  neglected.  Since  vm  -  v0  sin  0m,  for  m  -  1,2,  we 
substitute  values  of  sin  6m  directly  into  Eqn.  (13)  and  obtain  the  following 
approximate  relationship  between  the  strain,  e,  and  the  vibrational 
frequencies  vi  and  v2,  namely: 

e  -  (vi  -v2)/v*|  (14) 

We  note  that  compressive  in-plane  stress  is  usually  accompanied  by 
tensile  relaxation  in  the  out-of-plane  direction.  This  renders  an  optically 
isotopic  material  uniaxial.  In  the  geometry  we  are  using  for  the  ir  studies, 
the  electric  field  vector  is  always  in-plane,  and  that  the  ir  absorption 
spectrum  derives  from  modes  with  in-plane  polarizations.  As  such  the 
properties  of  these  modes  relfect  the  in-plane  stress,  and  as  shown  in 
Eqn.  (14)  also  yield  information  relative  to  the  in-plane  strain.  Since  the 
relationship  in  Eqn.  (14)  treats  differences  in  the  frequencies,  it  is 
important  to  determine  the  magnitude  of  the  terms  neglected  by  not 
treating  explicity  any  of  the  changes  in  the  Si-O  bond-length  [14].  These 
changes  come  from  differences  in  the  hybridization  of  the  oxygen  atom  2s 
and  2p  wavefunctions  as  20  is  changed.  As  20  increases,  the  amount  of  2s 
character  increases,  and  this  in  turn  means  that  the  Si-0  bond-length 
increases  as  the  Si-O-Si  bond-angle  is  decreased  [14].  Therefore  Eqn.  (14) 
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always  overestimates  the  strain.  We  now  demonstrate  this  in  Eqn.  (15), 
and  also  through  that  same  relationship  provide  a  quantitative  way  to 
determine  the  strain  to  the  next  level  of  approximation.  If  we  consider  the 
change  in  the  Si-0  bond-length  explicitly,  then  we  may  rewrite  Eqn.  (13) 
as: 


e  «  (2r0  sin0i  -  2[r0  +  5rj  sin  02)/2ro  sin  0*j,  (15) 

where  5r  is  the  change  in  the  Si-O  bond-length  associated  with  the  change 
in  the  Si-O-Si  bond-angle,  80  -  0-j  -02-  We  can  then  rewrite  Eqn.  (15)  in 
terms  of  vibrational  frequencies  and  obtain, 

e  -  (v-|-v2)/vi  -  (Sr/r0)(v2/vi).  (16) 

Combining  terms  defining  Sr  -  (Ar/A0)  80,  we  have, 

e  -  1  -  (1+[(Ar/A0)x(50/ro)])(v2/vi)).  (17) 

The  relationship  in  Eqn.  (17)  enables  us  to  estimate  the  strain  from  the  ir 
measurements,  provided  that  we  take  into  account  the  correction  term, 
[(Ar/A0)  80/ro](v2/vi),  which  we  can  obtain  from  the  results  presented  in 
Ref.  14.  We  determine  the  value  of  V2  from  ir  absorption  data  and  then 

determine  the  relative  strain  for  any  growth  temperature  and  oxide 
thickness  relative  to  the  reference  frequency  vi,  which  is  taken  to  be  the 

average  ir  frequency  for  a  thick  oxide  grown  at  1150°C.  The  bulk  material 
for  a  thick  oxide  grown  at  this  temperature  has  a  thermal  history  such 
that  the  growth  time  is  long  compared  to  the  viscoelastic  relaxation  time, 
and  the  region  in  which  there  is  a  stress  gradient  is  then  small  with 
respect  to  the  total  oxide  thickness.  This  is  displayed  in  Fig.  4  where  this 
behavior  is  clearly  evident  in  the  comparison  between  thermal  oxides 
grown  at  700°C  and  1000°C.  Note  again  the  growth  time  for  the  oxide 
formed  at  700° C  is  short  relative  to  t,  so  that  there  has  been  no  stress 
relaxation  in  that  oxide  film  during  its  growth  period. 

We  now  are  in  a  position  to  compare  stress  and  strain  in  the  SiC>2 

films,  and  thereby  deduce  an  elastic  constant  from  the  stress  and  ir  data. 
The  computation  of  a  renormalized  Young's  modulus,  E/(1-y)  [6],  for  Si02 

then  provides  a  check  for  the  validity  of  our  approach  regarding  the  way 
that  changes  in  the  ir  frequency  give  information  relative  to  strain  within 
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the  oxide  film;  y  is  Poisson's  ratio  for  Si02.  0.18  [17,18].  We  will  estimate 
E/(1-y)  from  the  data  for  the  oxide  film  grown  at  1150°C.  The  relationship 
between  Young's  modulus,  Poisson’s  ratio,  the  film  stress  and  the  relative 
strain  is  given  by: 

af  «  eE/(1  -y)  (18) 

where  cxf  is  determined  from  Eqn.  (1)  and  e,  the  relative  stress,  is  given  by 
Eqn.  (17).  <s\,  the  intrinsic  stress  at  the  Si/Si02  interface  is  determined 

from  the  intercept  in  Fig.  6  and  taken  to  be  4.6  x  109  dynes/cm2.  The  film 
stress,  of  is  obtained  by  combining  o;  with  the  thermally  generated  stress 

of  approximately  4  x  10^  dynes/cm2.  We  use  Eqn.  (17)  to  determine  the 
relative  strain.  We  take  0i  to  be  equal  to  150  degrees,  corresponding  to  a 
value  of  vi  -  1078.5  cm-1  (see  Eqn.  (7),  and  02  -  141.7  degrees, 
corresponding  to  v2  ■  1055  cm-1,  the  value  of  v,  extrapolated  from  the 
data  of  Fig.  4.  From  the  data  of  Ref.  14,  we  estimate  [(ArM0)/ro]  -  7.76  x 
10-4  degrees-1,  and  from  the  values  of  0-j  and  02  from  the  ir  data,  we  set 
50  «  8.3  degrees.  Proceeding  in  this  way  we  obtain  a  value  for  E/(1-y)  of 
approximately  (5.5  ±  0.5)  x  1011  dyne/cm2,  which  is  good  agreement  with 
the  quoted  value  of  approximately  8.5  xIO11  dyne/cm2  [17,18].  This 
comparison  demostrates  the  self-consistency  of  our  model,  and  the  use  of 
the  Si-Si  distance  as  a  measure  of  the  relative  strain. 

The  approach  used  in  the  determining  a  value  for  E/(1-y)  for  Si02 
from  the  total  stress  data  and  the  relative  strain,  as  calculated  from  Eqn. 
(17),  is  valid  for  the  film  grown  at  1000°C,  because  there  is  clearly 
viscoelastic  relaxation  of  the  interfacial  stress  in  the  growth  direction 
and  one  can  then  expect  linear  elastic  theory  to  apply.  This  in  turn 
justifies  the  use  of  the  renormalized  Young's  modulus  as  given  in  Eqn.  (18) 
[6].  On  the  other  hand,  the  time  scales  for  oxide  film  growth  and 
viscoelastic  relaxation  for  oxides  grown  at  700° C  are  such  that 
viscoelastic  relaxation  in  the  growth  direction  does  not  occur  during  film 
growth.  This  in  effect  means  that  the  700° C  data  can  not  be  treated  in  the 
same  way  using  a  renormalized  Young's  modulus  which  specifically  takes 
into  account  relaxation  in  the  z-direction.  To  illustrate  this  in  a  graphical 
way,  we  have  taken  intrinsic  stress  data  from  Fig.  6,  and  approximate 
strain  data,  estimated  from  Eqn.  (14)  and  Fig.  4,  and  plotted  these  in  Fig. 


9.  We  have  separately  analyzed  the  data  for  the  1000°C  and  the  700°C 
oxides,  and  the  results  of  these  linear  regression  analyses  are  included  in 
the  diagram.  As  expected  from  our  discussion  above,  we  find  that  the 
slope  of  the  1000°C  data  is  greater  than  that  of  the  700°C  data.  This  is 
interpreted  to  represent  the  differences  in  the  viscoelastic  relaxation 
phenomena  at  the  two  growth  temperatures;  i.e.,  the  slope  of  the  lOOOOC 
data  is  proportional  to  E/(1-y),  and  the  stress  extraplolates  to  zero  for 
zero  strain,  whereas  the  slope  of  the  700° C  data  is  greater,  and  the 
700°C  stress  does  not  extrapolate  to  zero  for  zero  strain.  Correcting  the 
1000°C  data  for  total  stress,  and  taking  into  account  the  bond-length 
changes  as  well,  brings  the  calculated  slope  into  good  agreement  with  the 
value  for  E/(1-y)  as  obtained  above.  On  the  other  hand,  the  regression 
analysis  fit  to  the  700°C  data  does  not  go  through  the  origin  indicating 
that  the  lack  of  viscoelastic  relaxation  in  the  z-direction  and  making  it 
not  possible  to  interpet  this  700°C  data  by  our  method  of  analysis.  The 
increased  slope  for  the  700° C  data  is  never  the  less  consistent  with  a 
compressive  in-plane  stress  and  a  lack  of  viscoelastic  relaxation  in  the  z- 
direction. 

V  SUMMARY 

By  way  of  summary,  we  make  the  following  points: 

(1)  The  compressive  stress  in  thermally  grown  Si02  films  is  a 

result  of  the  strain  induced  by  the  differences  in  the  molar  volumes  of  Si 
and  Si02  at  the  grown  interface.  This  leads  to  a  value  for  the  intrinsic 

growth  stress  at  that  growth  interface  that  is  independent  of  the  growth 
temperature  and  the  growth  rate.  Recent  experiments,  to  be  discussed  in  a 
later  publication  [20],  indicate  that  this  value  of  the  stress  at  the  Si/Si02 

interface  is  not  changed  during  high  temperature  [1 1 00° C]  thermal 
annealing.  This  means  that  the  effects  of  the  molar  volume  mismatch  at 
the  Si/Si02  interface  are  not  changed  via  annealing  even  at  temperatures 

where  substantial  viscoelastic  relaxation  can  occur  during  the  time  period 
of  the  annealing  cycle. 

(2)  The  extent  to  which  an  oxide  film  displays  a  stress  gradient 
depends  on  the  relative  time  scales  for  oxide  growth  and  viscoelastic 
relaxation.  Specifically,  for  oxides  grown  at  low  temperatues,  e.g.,  700°C 
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to  850° C,  the  viscoelastic  relaxation  times  are  very  long  with  respect  to 
the  growth  times,  and  stress  relief  in  the  growth  direction  via  a 
viscoelastic  mechanism  is  not  possible.  These  films  in  effect  have  a 
uniform  thermal  history.  This  is  characterized  experimentally  by  only 
small  differences  with  film  thickness  in:  (i)  the  intrinsic  stress;  and  (ii) 
the  ir  frequency  of  the  bond-stretching  vibration  which  provides  a 
measure  of  the  relative  strain.  On  the  other  hand,  for  oxides  grown  at,  or 
above  about  950° C  to  1000°C,  viscoelastic  relaxation  times  are 
comparable  or  shorter  than  the  growth  times,  and  substantial  viscoelastic 
relaxation  occurs.  This  leads  to  stress  gradients  and  the  associated 
variation  in  the  ir  vibrational  frequencies  which  reflect  an  accompanying 
variation  in  the  strain.  These  high  temperature  films  then  display 
evidence  for  a  non-uniform  thermal  history. 

(3)  The  viscoelastic  relaxation  times  estimated  from  stress  data, 
and  from  ir  data  are  approximately  the  same,  and  support  the  application 
of  a  Maxwell  model  for  viscoelastic  relaxation. 

(4)  The  computation  of  viscoelastic  relaxation  times  and  an  in-plane 
Young's  modulus  from  stress  and  ir  data  further  supports  the  use  of  the 
Si-Si  distance,  which  is  influenced  mostly  changes  in  the  Si-O-Si  bond- 
angle,  as  a  relative  strain  parameter. 

(5)  The  breakdown  of  a  universal  relationship  between  the  full- 
width  of  the  bond-stretching  ir  absorbance  and  the  ir  frequency  in  thin 
oxide  films  indicates  a  change  in  the  distribution  of  bond-angles  in  the 
immediate  vicinity  of  the  Si/Si02  interface.  The  result  described  above  is 

counter-intuitive  with  the  relative  distribution  of  angles  increasing  in 
the  vicinity  of  the  Si/Si02  interface  where  the  stress  and  strain  are 

pinned  by  differences  between  the  two  molar  volumes  at  the  growth 
interface. 

(6)  It  is  well-established  that  Si/Si02  interfaces  formed  at  lOOOQC 

and  incorporated  in  metal-oxide-semiconductor  (MOS)  structures  have 

relatively  defect  free  interfaces  with  the  densities  of  interfacial 
trapping  sites  as  deduced  from  C-V  measurements  generally  being  much 

less  than  about  5  x105 6 * * * 10/cm2  [i 9].  On  the  other  hand  Si/Si02  interfaces 

formed  at  700° C  generally  have  trapping  state  densities  that  are  more 
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than  an  order  of  magnitude  higher.  This  means  that  the  generation  of  the 
bonding  defects  that  lead  to  trapping  are  not  simply  related  to  either  the 
magnitude  of  the  intrinsic  growth  induced  stress  at  the  Si/Si02  interface, 

or  the  stress  gradient  in  the  immediate  vicinity  of  that  interface. 
Unfortunately,  our  study  does  not  supply  any  direct  information  relative  to 
other  possible  interfacial  defect  generation  mechanisms.  However,  our 
results  suggest  that  the  defect  generation  mechanism  may  in  part  related 
to  the  lack  of  stress  relief  in  the  low  temperature  thermal  oxides,  tox  < 

900°C. 

We  are  continuing  this  study,  and  have  set  up  experimental  studies 
to  examine  the  effects  of  annealing  on  thermally  grown  oxides  [20].  Our 
objective  is  to  vary  annealing  times  and  temperatures  for  oxides  grown  at 
low  temperature,  and  thereby  to  extend  the  observations  made  in  this 
paper  from  thermal  growth  to  thermal  annealing  as  well.  We  are 
particularly  interested  in  measuring  the  changes  in  the  stress  and  strain 
profiles  that  result  from  high  temperature  thermal  annealing,  and  in 
determining  the  extent  to  which  the  relative  times  scales  for  annealing 
and  for  viscoelastic  relaxation,  at  the  annealing  temperatures,  provide  a 
quantitative  basis  for  understanding  the  changes  brought  about  by  the 
annealing  process. 
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FIGURE  CAPTIONS 

Fig.1  IR  absorbance  for  Si02  bond-stretching  vibration.  The  figure  includes 
the  position  of  the  spectral  peak,  v,  the  full-width  at  half-maximum,  Av. 
and  an  insert  that  indicates  the  motion  of  the  oxygen  atom.  29  is  the  bond- 
angle  at  the  oxygen  atom  site. 

Fig.  2a  and  2b  v  (2a)  and  Av  (2b)  versus  oxidation  temperature,  respectively 
for  dry  oxides. 

Fig.  3  Av  versus  v  for  thermal  and  plasma  oxides,  including  both  types  of 
oxides  which  have  been  annealed.  The  line  is  a  regression  analysis  of  the 
data. 

Fig.  4  v  versus  oxide  thickness  for  thermal  oxides  grown  in  dry  oxygen  at 
700°C  and  1000°C.  The  solid  lines  represent  logarithmic  interpolations  of 
the  data.  The  approximate  intercept  for  both  sets  of  data  isl  055  cm-1. 

Fig.  5  Av  versus  oxide  thickness  for  thermal  oxides  grown  in  dry  oxygen  at 
700°C  and  1000OC. 

Fig.  6  Intrinsic  stress,  oj,  versus  oxide  thickness  for  thermal  oxides 
grown  in  dry  oxygen  at  700°C  and  1000°C.  The  extrapolated  stress  at  the 
Si/Si02  interface  is  estimated  to  be  about  4.6  xIO9  dynes/cm2. 

Fig.  7  Av  versus  v.  The  solid  line  is  the  regression  analysis  from  Fig.  3,  and 
the  points  are  for  thin  oxides  grown  in  dry  oxygen  at  700°C  and  1000°C. 
The  arrows  indicate  the  direction  of  increasing  film  thickness. 

Fig.  8  Av/v  versus  9  cot  9  [see  Eqn.  (10)].  The  solid  line  is  a  regression 
analysis  of  data  for  thick  oxides;  this  yields  a  slope,  A9/9  -  0.22.  The 
points  are  for  thin  oxides;  display  on  the  average  an  increased  slope  ■ 
0.47. 

Fig.  9  Intrinsic  stress,  o\,  versus  (vj-v2)/vi  as  determined  from  the  ir 
data  [see  Eqn.  (14)].  We  have  separately  analyzed  the  data  for  the  1000°C 
and  700°C  oxides,  and  the  lines  derived  from  a  linear  regression  analysis 
are  shown.  The  slope  for  the  1000°C  data  is  greater  than  that  for  the 
700°C  data.  In  order  to  obtain  estimates  of  the  elastic  constants  from 
these  slopes,  it  is  necessary  to  include:  (a)  the  thermal  contribution  to  the 
stress  data;  (b)  the  effects  of  viscoelastic  relaxation  in  the  z-direction; 
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and  (c)  the  contributions  of  bond-length  changes  as  a  function  bond-angle 
changes  to  the  strain  data. 
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Table  I  Estimated  Viscoelastic  Relaxation  Times,  x 

Oxidation  t,  Stress  data  x,  IR  data  x,  Viscosity 

Temperature  This  study  This  study  data,  Ref.  6 

700°c  ^>270  hours  332  hours  5,300  hours 


1000°C 


0.86  hours 


0.65  hours  0.20  hours 
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